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1  BACKGROUND 


In  the  study  of  the  Earth's  ionosphere,  there  are  many  problems  that 
require  a  detailed  knowledge  of  the  electron  distribution  function.  The 
electrons  are  produced  primarily  by  the  photoionization  of  t he  neutral  particles 
in  the  atmosphere  by  the  electromagnetic  radiation  from  tin  sun,  and  by 
secondary  impact  ionization  of  the  atmospheric  constituent  -  by  the  primary 
electrons.  The  electrons  undergo  various  collisional  interactions  with  the 
atmospheric  particles  and  a  steady-state  distribution  is  reached.  In  the  auroral 
atmosphere,  ionization  can  also  result  from  the  precipitation  of  the  high  energy 
(keV)  electrons  and  protons.  Once  the  electron  distribution  function  is  deter¬ 
mined,  ionospheric  properties  of  interest  can  then  be  calculated.  With  the 
advent  of  rocket-borne  and  satellite  electron  spectrometers  it  is  of  interest 
to  pursue  detailed  theoretical  studies  of  the  photoelectron  distribution  as 
well  as  of  the  electron  density  profile  in  the  auroral  ionosphere. 

The  Space  Data  Analysis  Laboratory  (SDAL)  of  Boston  College  was  contracted 
by  the  Ionospheric  Dynamics  Branch  (PHI)  of  the  Air  Force  Geophysics  Laboratory 
(AFGL)  to  develop  analytic  and  computer  techniques  for  theoretical  studies  of 
the  ionospheric  electron  distribution  function.  This  report  is  a  summary  of 
the  work  performed  under  the  auspices  of  SDAL  during  the  period  covered  by  the 


contract . 


PIIOTOELECTRON  I- LUX  IN  EARTH'S  IONOSPHERE  AT  ENERCl ES  NEAR  Till: 
PHOTO  IONIZATION  PEAKS 


In  this  work  Bolt zmann-Fokker-Planck  (BFP)  theory  is  used  to  analyze  the 
(isotropic)  photoelectron  flux  in  the  Earth's  ionosphere  at  energies  in  the 
vicinity  of  photoionization  peaks,  where  the  electron  impact  cross-sections 
are  slowly  varying  functions  of  energy. 

2.1  GREEN'S  FUNCTION  SOLUTION  FOR  THE  PHOTOELECTRON  FLUX 

In  the  steady-state,  local  approximation  an  equation  for  the  isotropic  part 
of  the  electron  distribution  function  is  obtained  by  integrating  the  BFP  equa¬ 
tion  over  the  angles  in  velocity  space.  The  equation  in  the  energy  representa¬ 


tion  is 


"£).*$  u  • 
pi  coll 


where  the  terms  on  the  right-hand  side  are  given  by  Eqs.  (7)  through  (15)  of 
Jaspcrse  (1976)  and  in  Jaspcrsc  (1975).  At  energies  in  the  vicinity  of  photo¬ 
ionization  peaks  the  momentum  transfer,  excitation,  and  ionization  cross- 
sections  are  slowly  varying  functions  of  energy  and,  as  a  result,  this  equation 
may  be  approximated  by 
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where  II  is  In  times  the  isotropic  part  of  the  electron  distribution  function 
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Here,  Q^.  are  the  momentum  transfer  cross  sections,  <}  j.  are  the  impact 

excitation  cross  sections  for  electronic  transitions,  0...  are  the  impact 

ionisation  cross  sections,  6.  =  2m/M  where  in  is  the  electron  mass  and  M.  is 

J  J  J 

the  neutral  particle  mass,  and  In  A  is  the  Coulomb  logarithm.  The  neutral 

particle  densities  are  denoted  by  n  with  an  appropriate  subscript  denoting  the 

type  and  state  of  the  neutral  particle,  the  neutral  particle  temperature  by  T 

(in  el') ,  and  ne  and  denote  the  analagous  electron  quantities.  The  angular 

brackets  denote  an  energy  average.  In  deriving  Eq.  (2)  rotational,  vibrational, 

fine  structure,  electron- i on ,  de-excitation  and  recombination  processes  were 

found  to  be  negligible.  The  production  of  secondary  electrons  by  electron 

impact  ionization  was  regarded  as  a  known  source  by  using  the  continuous  s lowing - 

down  approximation  and  was  included  in  S.  In  treating  impact  excitation  of  the 

electronic  states,  terms  of  the  form 
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where  are  the  threshold  energies,  were  also  neglected.  It  is  actually 
possible  to  solve  Eq.  (2)  with  these  terms  included  using  mathematical  methods 
more  general  than  the  ones  presented  in  this  paper.  These  terms  produce 
approximately  a  5%  effect  at  energies  in  the  vicinity  of  a  peak  produced  only 
by  photoionization.  These  more  general  results  will  be  published  elsewhere. 

Eq.  (2)  can  he  solved  if  we  are  able  to  find  the  Greta's  functions  which 
sat  isf ies 
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C (E , E ' )  =  6 (E-E ' ) 
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subject  to  appropriate  boundary  conditions.  Here,  A(E-E')  is  the  llirac  delta 
function.  The  photoelectron  flux,  4»(E),  is  then  given  bv 
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-1 


1/2 


/ 


dE'  G(E,E '  )  S(li') 


(4) 


I 


Using  standard  methods  we  obtain 
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where  the  boundary  conditions  that  G  >  0  as  E  -*■  00  and  that  G  have  no  rising 
solution  as  H  >-0,  have  been  imposed.  Here 


2  i.l  Two  rhotoeleetron  Sources 


The  first  photoelcctron  source  we  examine  is  that  due  to  a  continuously 
decreasing  source  and  an  isolated  unbroadened  photoionization  source. 

S(K)  =  So  exp(-l:/l-o)  +  S(  «CE-E1)  . 


(8) 


Using  hq.  (4),  the  photoelcctron  flux  is 
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Above  about  2(>0  km  electron-neutral  particle  collisions  are  d  .reusing  rapidly 

9 

and  electron-electron  collisions  dominate.  This  implies  that  b“>1ac  and  that 
(<j  1  /  1 0  '*  >  2c/i  n  .  Since  the  conti nous  source  decreases  slowly 

a |  '  1/b  ;  an  approximate  solution  is 
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At  high  altitudes  we  see  that  the  photoolcctron  flux  has  a  cusp  at  h  and  a  line 
shape  that  is  asymmetric  about  E  .  The  forward  and  backward  half  widths  at  the 
1/e  point  are 


FHW  =  T 

e 


(13) 


BHW  =  Y  n  /2c  ^  n  /n(0)  ,  (14) 

ee  c  e 

where  n(0)  is  the  atomic  oxygen  density.  At  low  altitudes  electron-neutral 

n 

particle  collisions  dominate  electron-electron  collisions  ;md  tT<4ac.  It 

1/7 

follows  that  a.~u0  =  (c/a)  £  a  *>  1/ii  and 

12  o 

(E)  ~  (4tt  )  ” 1  [l)/  exp(-E/Eo) 

+  -j  exp (-a | E-Ej |)  |  .  (15) 

The  photoelectron  flux  again  has  a  cost  at  E.  and  the  line  shape  is  symmetric 

1/a 

with  a  half  width  given  by  (a/c)  . 

The  second  nhotoelectron  source  we  examine  is  that  due  to  a  continuously 
decreasing  source  and  the  isolated  photoionization  source  that  is  broadened. 

S(F.)  =  Sq  cxp(-E/Eo)  +  S1  r (E)  ,  (10) 

where  r(E)  is  a  rectangular  function  of  height  1 / 2 r  and  of  width  21'  centered 
about  Ej.  The  photoelectron  flux  is  given  by  E<i .  (9)  with  fj  replaced  by  f  ’ 
where 
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Note  that  $  has  a  rounded  top  and  a  continuous  first  derivative  in  the 
vicinity  of  L(.  Note  also  it  is  asymmetric  at  high  altitudes  (a^o^)  and  sym¬ 
metric  at  low  altitudes  («  =m,,) .  At  altitudes  above  260  km  the  forward  half 
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for  simplicity  we  have  neglected  the  shift  in  the  location  of  the  peak.  For  T 
comparable  to  T  FllW  is  approximately  T  +F  -  T  In  2.  At  low  altitudes  the  half 
width  is 
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For  large  T(c/ a) 
(a/c)1/2  +  r. 


-omparid  to  1,  11W  is  approximately  given  by  (1-ln  2) 


2.2  Ki.sm/rs 

Using  the  model  ionosphere  same  as  that  given  in  .lasporse  (1077) f  and  using 
the  .Jacobin  (1977)  model  for  neutral  atmosphere,  we  have  calculated  the  theoreti¬ 
cal  photoelectron  flux  in  the  vicinity  of  the  27.2  eV  peak  at  three  altitudes. 
Results  of  these  c  ilculations  are  given  in  Jasperse  and  Smith  (1978).  The 

27.2  eV  peak  was  chosen  as  it  is  isolated  and  is  produced  only  by  photoionization 
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In  comparing  these  calculations  to  the  experimental  results  of  Doe  ring  et_  al^ 

( 1 97  <■> ) ,  we  see  that  the  observed  asymmetry  and  shift  of  the  photoioni  rat  ion 

peak  as  the  altitude  increases  is  qualitatively  explained  by  the  BFP  theory. 

REFERENCES 

Doering,  J.P.,  W.K.  Peterson,  C.O.  Bostrom,  and  T.A.  Potemra,  High  resolution 
daytime  photoelectron  energy  spectra  from  AE-E,  Gcophvs.  Res.  I.ctt.,  3^, 
129,  1976. 

Jacchia,  L.G.,  Thermospheric  temperature,  density,  and  composition:  New 
models,  Smithson.  Astrophys.  Observ.  Spec,  kept.,  375 ,  1977. 

Jasperse,  J.R.,  Electron  distribution  function  in  a  nonuniform  magnetized 
weakly  photoionized  gas  application  to  a  model  ionosphere,  AI-'CRL  Rep. 
TR-75-0266 ,  1975. 

Jasperse,  J.R.,  Boltzmann-Fokker-Planek  model  for  the  electron  distribution 
function  in  the  Earth's  ionosphere,  Planet.  Space  Sci . ,  24 ,  33,  1976. 

Jasperse,  J.R.,  Electron  distribution  function  and  ion  concentrations  in  the 
Earth's  lower  ionosphere  from  Boltzmann-Fokker-I’lanck  theory.  Planet . 
Space  Sci . ,  25,  743,  1977. 

Jasperse,  J.R.  and  Ed  R.  Smith,  Geophys.  Res.  Lett.,  5_,  843,  1978. 


7 


3  COLLISION  OPERATOR  STUDY 

We  have  completed  ;i  systematic  study  ol'  collisional  operators  based  on 
the  "velocity-ratio"  expansion  technique  introduced  originally  by  Scluillcr  and 
Wilhelm1.  This  technique  has  been  developed  extensively  by  Sdgur  and 
Lerouvi llois-Gai Hard2 ,  and  is  applicable  to  all  types  of  collisional  processes 
involving  light  and  heavy  particles.  It  is  particularly  convenient  when  con¬ 
sidering  the  impact  of  electrons  with  atoms,  molecules,  or  ions.  This  calcu- 
lational  procedure  considers  systematically  the  dominant  effects  of  the  number 
density,  temperature,  drift,  velocity-correlation  tensor  (or  stress  tensor)  of 
the  heavier  species  on  the  collisional  operators  of  the  transport  equation  for 
the  distribution  function  of  the  lighter  species. 

In  our  study,  we  have  considered  the  effects  of  ionization,  dissociation, 
excitation,  de  excitation,  recombination,  and  Coulomb  forces  on  the  collisional 
operators.  the  following  is  a  summary  of  the  specific  results  of  our  investi- 
<>at  i  ofi . 

3.1  COLLISION  OPERATORS  VIA  SCATTERING  PROBABILITIES 

In  considering  collisions  between  electrons  and  neutral  species,  we  found 
it  convenient  to  unploy  the  "scattering  probability"  formulation  of  the 
Bol t zmann- 1 i ke  operator.  This  idea  was  discussed  by  Wnldmann3  and  used 
extensively  bv  He;  logo1*  and  Scgur  and  Lcrouvi 1 loi s-Ga i 1  lard2 .  It  allows  us 
to  treat  all  Bo  1 1 .  inann- 1 i ke  collision  operators  from  a  single  unified  viewpoint. 
We  have  used  this  procedure  to  consider  the  inelastic  scattering  between 
electrons  and  neutral  species  and  electron-neutral  impact  ionization  processes. 


3.1.1  Inelastic  Scattering  Between  Electrons  and  Neutral  Species 

We  have  considered  the  inelastic  scattering  of  neutral  particles  (atoms) 
and  molecules)  with  colliding  electrons.  The  atoms  or  molecules  are  allowed 
to  undergo  excitation  and  tie-excitation  processes.  Principle  of  detailed 
balance  is  used  to  obtain  the  final  form  of  the  collision  operator.  The  dis¬ 
tribution  functions  of  the  neutral  particles  are  assumed  to  he  given  but  their 
functional  forms  an  he  quite  arbitrary. 


i)  We  have  derived  the  eollisional  operator  for  excitation  processes 
to  second-order  (in  velocity-ratio)  for  ;in  isotropic  neutral 
distribution  function.  This  result  is  the  j;cncr;i  1  i  zat  i  on  of  the 
classical  expression  developed  by  Holstein’.  In  the  original 
expression  of  Holstein,  the  collision  operator  depends  only  on  the 
number  density  of  the  colliding  neutral  species  and  the  ili fferent ial 
cross  section.  In  our  second-order  calculation,  the  neutral  tempera¬ 
ture  is  also  included  in  the  final  expression. 

We  have  written  the  operator  in  an  expanded  form  by  writing  the 
electron  distribution  function  in  spherical  harmonics  of  the  angles 
of  the  (electron)  velocity  vector.  The  first  two  moments  of  the 
operator  have  been  obtained  explicitly.  They  correspond  to  the 
isotropic  part  and  the  first  anisotropic  contribution  of  the  colli¬ 
sion  operator.  As  a  check,  we  have  shown  that  in  the  limit  of  zero 
excitation  energy,  the  isotropic  part  of  the  collision  operator 
reduces  to  the  classical  Fokker-Planck  form  of  the  elastic  collision 
operator.6  Thus,  the  elastic  eollisional  contribution  is  a  second- 
order  effect  in  terms  of  the  ratio  of  the  velocities  of  the  colliding 
neutral  particle  and  electron.  It  is  not  contained  in  the  original 
Holstein  expression. 

ii )  We  have  derived  the  corresponding  collision  operator  to  second 

order  for  de-excitation  processes  for  an  isotropic  neutral  distribu¬ 
tion  function.  This  is  accomplished  by  applying  the  "principle  of 
detailed  balance"  to  the  excitation  and  de-excitation  processes. 

By  expanding  both  the  excitation  and  de-excitation  operators  in 
Taylor  series  in  terms  of  the  excitation  energy  and  summing  the 
resulting  expressions,  we  obtain  a  generalized  version  of  the 
Ginsburg-Gurevich-Al 1 i s7  form  of  the  isotropic  part  of  the  collision 
operator.  In  this  calculation,  we  have  assumed  quasi -equi librium 
of  the  neutral  species  and  applied  the  Saha  equation.  Thus  the 
effect  of  the  neutral  temperature  is  introduced  in  the  collision 
operator  through  both  the  second-order  velocity-correlation  func¬ 
tions  of  the  neutral  species  via  the  velocity-ratio  expansion 
technique  and  also  through  the  expansion  of  the  Saha  expression. 

Our  expression  reduces  to  the  classical  Fokker-Planck  form  if  only 
the  terms  to  zeroth-order  in  velocity-ratio  and  second-order  in 
excitation-energy  arc  retained.  [We  note  that  this  expression  still 
contains  the  neutral  temperature  because  of  the  Saha  equation).  To 
second  order  in  velocity-ratio,  however,  the  structure  of  this 
collision  operator  is  more  complicated  than  the  original  Fokker- 
Planck  form.  We  must  be  careful  of  not  double- summing  the  elastic- 
terms.  Now  the  operator  contains  both  first  order  and  higher- 
order  terms  in  excitation  energy,  as  well  as  third  and  fourth  order 
derivatives  of  the  electron  distribution  function. 

i i i )  We  have  derived  the  collision  operator  to  first  order  (in  velocity- 
ratio)  for  inelastic  impact  for  an  anisotropic  neutral  distribution 
function.  This  expression  contains  both  the  number-density  and 
dr ift -velocity  of  the  neutral  species.  We  have  expanded  this 
operator  in  spherical  harmonics  of  the  angles  of  the  (electron) 
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velocity  vector.  Since  t  lie  spheric;)  I  harmonics  ;iro  no  longer' 
eigenfunctions  of  the  collision  operator,  we  f  ind  that  the  isotropic 
portion  of  the  collision  operator  depends  on  the  anisotropic  part 
of  the  electron  distribution  function  and  vice  versa. 

3.1.2  Electron-Neutral  Impact  I  ionization 

We  have  considered  the  collision  process  of  electrons  and  neutrals  such 
that  the  neutrals  are  singly  ionized  by  emitting  secondary  electrons.  For 
simplicity,  we  have  assumed  that  the  neutral  distribution  function  is  isotropic. 
Thus,  we  shall  neglect  the  effect  of  neutral  drift  on  impact  ionization. 

The  collision  operator  is  found  to  depend  on  a  double  differential  cross 
section  as  discussed,  c.g.,  in  Mott  and  Massey8.  To  second  order  (in  velocity- 
ratio),  our  expression  depends  on  the  number  density,  and  temperature  of  the 
neutral  species.  Vc  have  expanded  the  operator  in  spherical  harmonics  of  the 
neutral  species.  We  have  expanded  the  operator  in  spherical  harmonics  of  the 
(electron)  velocit  ■  vector.  The  resulting  expressions  depend  on  the  partial 
range  integrals  ol  the  double  differential  cross  section.  Because  of  the 
impact  and  secondary  electrons  are  indistinguishable  after  the  collisional 
process,  a  factor  >f  two  (2)  is  introduced  in  some  of  the  partial -range 
integra  Is . 

To  zeroth  order  (in  velocity-ratio),  the  isotropic  portion  and  t he  first 
anisotropic  moment  of  this  collision  operator  reduce  to  the  expressions 
derived  previously  by  Jasperse.8 

3.1.3  Collision  Operator  for  Dissociative  Recombination 

Wo  arc  primarily  interested  in  the  rccombinat ive  processes  of  electrons 
with  molecular  ions  such  that  the  molecules  are  dissociated  into  two  atoms 
(or  molecules)  upon  collision.  The  inverse  of  such  a  process  does  not  involve 
collisions  with  the  electrons  and  tiiis  portion  of  the  collision  operator  is 
independent  of  th<  electron  distribution  function.  It  is  essentially  a  "source 
term".  Thus,  it  will  be  convenient  to  consider  the  collision  operator  in  two 
distinct  parts;  me  involving  the  forward  process  which  is  a  linear  function 
of  the  electron  distribution  function  and  t he  other  a  source  term. 

'io  consider  the  part  of  the  collision  operator  involving  the  forward 
rccombinat i vc  proi css  of  electrons  with  molecular  ions,  we  expand  the  differen¬ 
tial  cross  section  in  terms  of  the  "ion-electron  velocity-ratio".  The  result 
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is  a  simple  expansion  of  the  collision  operator  for  rccombi native  processes. 
The  effects  of  the  number  density,  drift,  and  temperature  <  f  the  ion  species 
come  in  naturally.  To  the  zeroth  order  of  this  expansion,  our  result  agrees 
with  the  classical  expression.  (Sec  e.g.,  Jasperse9.) 

3.1.4  Effects  of  Ion  Drift  on  the  Fokker-Planck  Operator  'or  Electron  - 
Ion  Collisions 

It  is  known  that  the  Fokker-Planck  form  of  the  collision  operator  for 
charged  particles  can  be  deduced  from  the  Boltzmann  equation10,  Liouville 
Theorem  or  Markov  processes12.  In  considering  electron- ion  collisions  and  the 
collision  operator,  it  has  been  customary  to  assume  a  Maxwellian  distribution 
for  the  ions.  It  is  quite  straightforward  to  generalize  the  expressions  to 
arbitrary  ion  distributions  using  the  velocity-ratio  expan  ion  technique.  We 
have  systematically  derived  the  expanded  operator  to  include  the  effects  of 
the  ion-drift  and  ion-stress-tensor .  We  have  again  splitt>-d  the  resulting 
operator  into  an  isotropic  part  and  the  first  order  anisotropic  moment.  Our 
results  reduce  to  the  classical  expression  in  the  limit  of  zero  ion  drift  and 
isotropic  ion  distribution. 

3.2  INTEGRATION  OF  THE  ELECTRON  TRANSPORT  EQUATION  IN  APPLIED  ELECTRIC 
AND  MAGNETIC  FIELDS 

We  have  extended  the  previous  calculations9'1  of  the  electron  distribution 
function  in  the  collisional  region  of  the  ionosphere  to  include  the  effects  of 
the  local  electric  and  magnetic  fields.  In  a  local  approximation,  we  shall 
assume  that  the  E,B-fields  are  given.  We  found  that  the  integration  of  the 
electron  transport  equation  can  be  accomplished  conveniently  if  we  project 
the  electron  distribution  function  in  three  (3)  special  directions: 

i)  the  direction  of  the  electric  field, 

ii)  the  direction  of  the  magnetic  field,  and 

iii)  a  direction  normal  to  both  the  electric  and  magnetic  field. 

We  found  that  if,  for  the  isotropic  portion  of  the  c<  llision  operator,  we 

i)  neglect  recombination  and  attachment, 

ii)  use  the  Ginsburg-Gurcv ich-Al i is  form  of  the  inelastic  collision 
operator,  or  expand  the  inelastic  operator  |to  zero  order  of  the 
velocity  ratio]  in  Taylor  scries  in  terms  of  the  excitation  energy, 
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iii)  linearize  the  electron-electron  collision  term, 
ivj  neglect  l  lie  neutral  and  ion  drifts, 

v)  approximate  electron-neutral  impact  ionization  as  an  excitation 
process  and  expand  in  Taylor  series  in  terms  of  the  excitation 
energy, 

vi )  consider  photo  ion i zat ion ,  and  the  secondary  electron  production 

due  to  electron-neutral  impact  (estimated  by  the  continuous  slowing 
down  approxiroat i on)  as  known  sources,  and  use  a  combined  collision 
frequency  for  the  anisotropic  part  of  the  collision  operator,  the 
electron  transport  equation  can  be  integrated  in  terms  of  simple 
quadratures.  Interestingly,  these  quadratures  can  be  integrated 
explicitly  for  simple  but  realistic  variations  of  the  cross 
sections  ns  a  function  of  the  electron  energy. 

We  note  that  the  above  integration  is  accompli  shed  by  simultaneously 
considering  the  isotropic  and  anisotropic  portions  of  the  transport  equation 
for  the  electrons. 

3  .5  INTEGRATION  OF  THE  ELECTRON  TRANSPORT  EQUATION  IN  THE  ABSENCE  OF  AN 

ELECTRIC  FIELD 

Ln  the  absence  of  the  elhctric  field,  we  can  integrate,  in  local  approxi¬ 
mation,  the  electron  transport  equation  perturbatively .  If  we  use  the  same 
type  of  approximations  suggested  in  Section  2  for  the  collision  operators, 
the  isotropic  part  of  the  electron  transport  equation  can  be  integrated 
explicitly.  This  expression  is  probably  only  valid  in  the  electron  energy 
range  JU  eV,  but  should  he  quite  accurate  for  the  pbotoioni zat ion  peaks  regime 
and  higher  electron  energy  ranges. 

The  anisotropic  part  of  the  electron  transport  equation  in  terms  of  a 
single  combined  collision  frequency  can  also  be  integrated  in  the  presence  of 
an  applied  magnetic  field. 

In  an  earlier  study,  .laspcrso  and  Smithgc  calculated  the  photoelectron 
flux  in  the  vicinity  of  photo  ion i zat ion  peaks  in  terms  of  a  local  Green's 
function  solution  for  the  isotropic  part  of  the  electron  distribution  function. 
We  have.-  applied  the  same  idea  to  calculate  the  anisotropic  effects  in  the 
presence  of  an  applied  magnetic  field  near  photoionization  peaks. 

3.4  SUMMARY 

We  have  compiled  a  systematic  study  of  the  collision  operators  that  appear 
in  the  electron  transport  equation  using  the  velocity-ratio  expansion  technique. 
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Processes  involving  inelastic  impact,  excitation  and  dc-ex<  itat ion,  ionization, 
dissociating,  recombination  and  Coulomb  collisions  have  be.  n  considered.  We 
believe  this  research  is  the  most  comprehensive  and  systematic  investigation 
of  the  collision  operators  for  the  scattering  of  light  particles  by  heavy 
particles  at  the  present  time. 

Upon  making  some  plausible  approximations  we  are  able  to  integrate  the 
electron  transport  equation  in  the  presence  of  electric  and  magnetic  fields  in 
the  local  approximation. 

We  are  in  the  process  of  preparing  two  separate  manuscripts  reporting 
these  research  findings. 
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4  I: l J iCTROSTAT I C  PLASMA  INSTABILITIES  01-'  li -REGION  I’JIOTOELHCTRON 
1)1  S  I  R  I  BU  T  IONS 

4  . 1  INTRODUCTION 

An  ionis'd  region  exists  in  the  Earth's  upper  atmosphere ,  which  extends 
from  about  sixty  kilometers  above  the  Earth's  surface  to  several  thousand 
kilometers  and  beyond.  At  low  latitudes  during  day-light  hours,  this  ionized 
region  results  primarily  from  the  photoionization  of  the  neutral  atmospheric 
constituents  by  eli ctromagnetic  radiation  from  the  sun.  In  the  study  of  this 
ionized  region  there  are  many  problems  of  interest  which  require  a  detailed 
knowledge  of  the  electron  energy  distribution.  Once  the  distribution  function 
is  known,  various  important  physical  quantities,  such  as  the  volume  emission 
rates  for  the  air  glow  radiation,  electron  heating  rates,  and  electron 
temperature  profiles,  can  be  calculated. 

Theoretically,  the  electron  distribution  function  can  be  determined  by 
invoking  appropriate  energy  balance  between  the  photoelectron  production  and 
the  losses  due  to  various  inelastic  processes.  Several  authors  have  investi¬ 
gated  this  problem.  Most  recently,  .Jasperse1  has  developed  a  complete  kinetic 
theory,  using  the  Boltzmann-Pokkcr-Planek  method,  in  the  sense  that  it  deter¬ 
mines  the  electron  energy  distribution  and  the  ion  concentrations  in  the  lower 
ionosphere  (E  region)  scl f-consi stently ,  once  the  boundary  conditions  and  the 
model  ionosphere  are  specified.  The  low  energy  spectrum  of  the  isotropic  part 
of  the  distribution  function  is  Maxwellian  with  temperatures  a  few  hundredth 
of  an  eV.  In  the  upratherma 1  part  of  the  spectrum,  most  interesting  feature 
occurs  between  2  and  4  eV  energy.  In  this  energy  range,  the  distribution 
function  has  a  minimum  at  about  2.3  eV  (Pig.  I).  This  is  explained  theoreti¬ 
cally,  since  the  c  oss  section  for  the  electron- impact  excitation  of  vibration¬ 
al  states  of  N,  ha;  a  maximum  at  2.3  eV.  Beyond  2.3  eV,  the  photoelectron 
flux  rises  sharply  to  a  maximum  at  about  4  eV,  as  the  excitation  cross  section 
decreases  sharply.  Beyond  this  energy  the  flux  falls  off  as  the  cross  section 
for  the  excitation  of  metastable  states  of  atomic  oxygen  increases.  The 
theoretical  calcul  itions  (Pig.  1)  show  that  this  minimum  in  the  spectrum  is 
more  and  more  pron  mneed  at  or  below  altitudes  of  130  km.  Above  130  km,  the 
valley  starts  to  I"-  filled  up,  the  peak-to-val ley  ratio  decreases,  until  at 
or  above  210  km  tie-  structure  disappears  completely.  This  disappearance  can 
be  attributed  to  the  depletion  of  N,  as  well  as  to  the  gradual  smoothing  process 
arising  from  electron-electron  collisions  as  the  altitude  increases. 
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Measurements  of  the  electron  energy  distribution  bv  Doerinj  £t_  and,  most 
recently,  by  McMahon  and  Heroux3,  who  studied  specifically  t lie  2-5  eV  energy 
range  with  improved  energy  resolution  of  the  apparatus,  art  in  good  agreement 
with  the  theoretical  calculations  of  Jaspersc1  at  and  above  170  km.  Below 
170  km,  the  calculated  values  of  the  peak-to-val ley  ratios  by  Jaspersc  is 
somewhat  larger  than  the  observed  values.  The  observation;  are  in  most  striking 
disagreement  with  the  theory  below  130  km,  where  they  show  plateaus  in  the 
distribution  functions  in  the  2-4  eV  energy  range. 

This  discrepancy  between  theory  and  observations  suggests  that  in  the  low 
altitude  regions  (100-170  km)  collisional  processes  along  cannot  account  for 
the  observed  photoelectron  distributions.  It  is  well-known  that  a  homogeneous 
plasma  in  a  magnetic  field  with  isotropic  distribution  functions,  can  he 
unstable  if  a  population  of  high  energy  particles  is  also  present.  This  is 
precisely  the  situation  in  the  lower  ionosphere,  and  it  is  expected  that  the 
plasma  instability  will  produce  the  anomalous  diffusion  in  the  velocity  space 
through  wave-particle  interaction,  which  in  turn  will  flatten  the  distribution 
functions  in  the  2-4  eV  energy  range.  With  this  in  mind,  \.e  consider  here 
excitation  of  electrostatic  instabilities  in  the  ionosphere  collisional  plasma 
by  the  suprathermal  electrons  near  the  4-eV  maximum  (fig.  !).  In  particular, 
we  consider  the  two  most  important  instabilities,  namely,  t he  upper  hybrid 
instability  and  the  electron  cyclotron  instability;  and  discuss  the  linear 
growth  rates  and  the  wave  spectra  as  functions  of  altitude  Bloomberg1"  bus 
investigated  the  upper  hybrid  instability  using  the  theorei  ical  electron  dis¬ 
tribution  function  of  Oalgarno  et  a  1 b .  His  calculations  show  that  the  instability 
is  inoperative  at  or  below  loO  km  due  to  large  collisional  damping  rate.  Since 
Jasperse1  gives  more  accurate3  theoretical  electron  distribution  functions  and 
since  we  now  have  more  accurate  data  for  determining  the  relevant  collision 
frequencies,  it  is  worthwhile  to  reexamine  the  upper  hybrid  isntnbility. 

In  Sec.  4.2,  we  present  the  mathematical  formulation  leading  to  the  dis¬ 
persion  relation.  In  Sec.  4.3,  we  solve  the  dispersion  relation  and  determine 
the  various  growth  rates,  and  the  spectrum  of  the  unstable  waves. 

4.2  MATHEMATICAL  FORMULATION 

Since  we  intend  to  study  the  instabilities  driven  by  the  suprathermal 
electrons  near  4-eV  maximum  (Fig.  1),  we  assume  that  the  equilibrium  electron 
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popu  1  ;i t  i on  is  composed  of  only  two  parts:  ;i  Maxwellian  (thermal)  population, 

and  a  siipratlicrinal  oik-  with  monotonical  Jy  decreasing  energy  profile  peaked  at 

4  eV.  This  is  indicated  by  dotted  lines  in  l-'igure  1.  In  doing  so,  we  have 

ignored  the  high  energy  photoi oni zat ion  peaks  in  the  tail  of  the  actual  electron 

distribution  function.  We  further  assume  that  the  positive  ions  form  a 

neutralizing  background  of  constant  density.  This  is  justified  since  we 

consider  waves  with  frequencies  u>w  .  where  u>  is  the  electron  plasma  fre- 

1  pe  pe 

quency.  Then  the  small  amplitude  perturbation  of  the  electrons  is  governed  by 
the  following  equations: 


3f . 

■L  +  v*Vf. 
4t  -  ~  J 

c 

me 

(V 

V  *  11  -  -4  FTC  k 

fi\ 

d3V 

_ L  _  —  p 

3v  m  ~ 


C(f.) 


where  B  is  the  earth's  magnetic  field,  which  is  assumed  to  be  uniform, 

f.(r,v,t)  is  the  perturbed  electron  distribution  function,  1  . (vj  is  the 
l  -  -  1  oj  ' 

equilibrium  election  distribution  function  isotropic  in  velocity  space,  and  1: 
is  the  electric  field  associated  with  the  density  perturbat ion.  The  index  j 
labels  the  two  kinds  of  electron  population. 

Hie  collision  term,  l',(f.J,  in  kq.  (1)  represents  the  electron-electron, 
electron- ion  and  the  elect ron -neutra 1  collisions.  Explicit  forms  of  these 
collision  terms  can  be  found  in  Kef.  (>.  in  a  weakly  ionized  gas  such  as  the 
low  altitude  ionosphere  of  our  interest,  the  frequencies  of  electron-electron 
and  electron-ion  collisions  are  small  in  comparison  with  the  neutral  and  wc 
need  consider  on  1 \  the  electron-neutral  collisions.  Furthermore,  inelastic 
electron-neutral  >.  >llisions  are  much  less  frequent  than  the  elastic  ones. 
However,  t  lie  complexity  of  the  exact  elastic  collision  term  makes  the  analysis 
intractable.  To  avoid  this,  we  use  a  B-(i-k  type7 
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is  the  perturbed  electron  density,  n  is  the  equilibrium  density,  and  v.  is 
the  sum  of  the  frequencies  ol  elastic  collisions  between  the  j  type  electrons 
and  the  various  neutral  particles  (N.,,  O,,  (),  etc.)  present  in  t  he  ionosphere. 


)(> 


This  particular  choice  of  the  collision  term,  which  conserves  only  the  number 
of  particles,  is  adequate  to  describe  the  electron-neutral  collisions. 


We  notice  that  Eq.  (1)  can  be  written  as 


df . 
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dt  j  j  m 
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where  d/dt  is  the  convective  derivative  taken  along  the  unperturbed  electron 
orbit  in  the  magnetic  field  B  Integrating  Eq.  (4)  with  respect  to  t  we 


obtain 
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F  .  (v ') 
oj 


(5) 


where  jr'(t')  and  v^'(t')  are  the  unperturbed  electron  orbits.  Since  the 
collision  frequency  is  much  less  than  the  cyclotron  frequency,  the  electron 
orbits  in  lowest  order  are  given  by 


x'-x 


s  i  n4> 


s  in  (4>-U  t) 


J_  I 

y'-y  =  -jj  cos*  -  jp  cos (*-Slt ) 

Z'-Z  =  -  VUT 

VX  =  vi-  C0S^"fiT^ 

Vy  =  vx  sin(4>-fiT) 

v '  =  v,, 
z  t\ 

where  t  =  t-t',  and  we  have  chosen  a  cylindrical  coordinate  system  in  the 
velocity  space  with  the  z-axis  along _Bq,  so  that  v  =  (v^cos$,vAsin$ ,v,%) .  The 
orbits  are  taken  to  satisfy  r'(t'=t)  =  r,v'(t'=t)  =  v,  and  SI  =  eB  /me  is  the 

~  >  -S  o 

electron  cyclotron  frequency. 

We  use  ^E  =  -V*,  where  *  is  the  electrostatic  potential,  and  choose 
perturbations  of  the  form  f  =  f  exp(ik*v- iwt) .  Then 
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4.3  STABILITY  AN. LYSIS 

4.3.1  Uesonant  I  pc  I  nst  al>  i  I  i  t  v 

Here  we  look  for  instabilities  that  arise  due  to  wave-particle  resonant 
interactions.  Th<  waves  are  •aipport oil  l>v  the  Maxwellian  electrons,  while  they 


IS 


arc  driven  unstable  by  the  suprathe rma I  electrons  near  tin  I  eV  peak,  which 

are  in  resonance  with  the  waves.  The  resonance  condition  i  s  to  -  k,.v,  =  nil.  to 

r  't  h  r 

being  the  real  part  of  the  frequency.  We  first  consider  the  properties  of  the 
waves  by  solving  the  dispersion  relation  without  the  suprathermal  electron 
term,  i.e.,  we  solve  the  equation 
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4tte 


2  D 
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where  the  subscript  m  refers  to  Maxwellian  electrons,  for  Maxwellian  electrons 
with  density  n^  and  temperature  T  ,  the  integrations  in  N  and  D  can  be  carried 
out  in  a  straightforward  manner  and  the  results  are 
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Here  Z(p^)  is  the  plasma  dispersion  function,  ln  is  the  modified  Bessel  function 
of  first  kind,  =  (w+iv^-nil)  /  [  |  k  |  (2T^/m) 1  ^ ,  b  =  k^  T^/ .  We  shall 
consider  waves  for  which|p  |>>1  for  all  n.  This  enal>lcs  us  to  ignore  the 
Landau  damping.  Using  the  asymptotic  expansion  for  the  Z  function0  and  keeping 
the  significant  leading  terms,  we  get 
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We  consider  waves  with  k()/k^<<l,  and  since  v 
of  Hqs.  (14)  and  (15)  , 
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2  1/2 

where  w  =  (4nn  e  /m)  is  the  electron  plasma  frequency,  lixpanding  in  v  /w<<1 
pe  o  m 

and  keeping  only  the  leading  terms,  we  can  write  bq.  (16)  in  the  form 
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where 
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If  we  write  uj  =  w  *  iui.  ,  where  i».  w  in  consistency  with  v  <<  m,  then  u  , 
r  i  i  r  1  m  r 

is  determined  from 


>  (k  ,k  ,u  )  =  0 
and  w.  is  given  by 
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In  deriving  liq.  (2ft),  which  gives  the  col lisional  damping  of  the  modes,  we 
have  made  use  of  the  dispersion  relation  (19). 


The  roots  of  the  dispersion  relation  in  the  limiting  case  k(l=0  was  first 


discussed  by  Bernstein.10  They  depend  on  the  values  of  the  parameters:  b  and 
2  2 

ft  /cu  t.  In  the  cold  plasma  limit  h< 1 ,  and,  furthermore,  over  the  ionospheric 
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region  of  our  interest  (100  170  km)  ,  ft“/u  •'<1.  Then,  the  lowest  frequency 
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Ihe  next  higher  fnquency  mode  is  the  so-called  upper  hybrid  mode: 
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(22) 
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The  higher  frequency  modes  (just  above  or  just  he  low  the  cyclotron  harmonic 

9  ? 

depending  on  the  value  of  u  /S2“)  are  given  by  (n-3,4 . i 

o»n  =  nfl  [1  -  (Kb11'1)]  ,  -  (n2- 1 )  Si2 


-  nfi  [1  +  0 (bn~ 1 )  ]  ,  (o2  <  (n2-l)  U2  ( 

pc 

These  so-called  Bernstein  modes10  are  separated  by  frequent y  intervals  where 

2 

there  is  no  propagation,  i.e.,  kx  *-  0.  In  Eqs.  (21)  -  (23)  there  are  small 

9  9 

corrections  of  the  order  of  k"/k^  <<  1,  which  have  been  omitted  for  the  sake 
of  convenience.  The  corresponding  damping  rates  are  found  to  he 
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respectively  for  the  modes  uk,  co,,,  and  oj  . 

2  H  n 

In  order  to  find  the  destablizing  effect  of  the  suprathermal  electrons 
(denoted  by  a  subscript  h) ,  we  now  consider  the  complete  dispersion  relation: 

[L.H.S.  of  Eq.  (17)]  +  (4-ne/k2)  (Nh/D  )  =  0  (1 

where  and  are  to  be  calculated  from  Eqs.  (8)  and  (9)  respectively.  We 
write  w  =  +  iw.  +  iy,  where  ok  is  the  collisional  damping  rate  given  by 

Eq.  (20),  and  y  is  the  growth  rate.  Since  dv  =  n^  '  nQ,  we  can  assume 

the  the  contribution  of  the  suprathermal  electrons  to  the  dielectric  response 
function  and,  hence,  the  growth  rate  is  small.  Then,  as  before,  we  obtain 


C(kX,kU’Wr'i  =  0  ’ 


which  del  ermines 
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where  c  (kj_ ,k(l,u>)  is  the  same  as  in  Eq.  (18).  In  doing  so  we  have  neglected 

the  small  contribution  to  m  from  Re(N,/D,  ). 

r  h  h 

We  first  recognize  that  the  collisions  of  the  suprathermal  electrons  with 
the  neutral  species  will  somewhat  detune  their  resonance  with  the  waves.  As 
a  result,  y  will  be  less  than  its  eol I isionless  value.  However,  the  reduction 
is  found  to  be  negligibly  small,  c.g.,  y  =  y  [1  -  OCv^/Q)],  where  <<  1, 

Y  being  the  collisionless  growth  rate.  This  can  he  easily  verified  by 
carrying  out  the  calculation  with  a  simple  functional  form  for  F  ^(v),  e.g., 

=  (n^^/^v”)  «lv-v  ),  being  the  speed  corresponding  to  the  peak  energy. 
So,  the  dominant  damping  is  that  due  to  the  collisions  of  the  Maxwellian 
electrons  and  is  given  by  Eqs.  (24)  -  (2(>) .  For  the  sake  of  convenience,  we 
shall  omit  from  and  Then  1>^  becomes  unity,  and  since  F  ^  is  isotropic 

in  velocity  space  N  becomes 
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Substituting  this  into  Eq.  (28)  and  using  the  formula 
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where  the  symbol  P  denotes  the  Cauchy  principal  value,  we  finally  get 
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Wc  change  over  lo  the  energy  represent  at  ion  using,  G^^U-J  -  ( I;  "/m  ) 

Foh(v),  so  that  the  normalization  condition  if  4n  J“„h<i:>  Jl  */'o h‘vl  J-v  ■  '• 
Here  E  =  mv  /2  is  the  energy  variable.  In  the  energy  representation,  Eq.  (30) 
becomes 
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where  Eq  =  mv^/2  is  the  peak  energy  of  the  suprathermal  elvetrons,  K(^  k^v^/U, 
and  E^  "  [  (wr/fl-n)/K4(]  ^ .  Integrating  by  paits  and  using  a  Bessel 

function  identity  we  rewrite  Eq.  (31)  in  the  more  convenient  form 
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The  behavior  of  the  functions  in  the  integrand  of  Eq.  (32)  is  shown  in  Figure 
2.  It  is  clear  that  to  obtain  large  contribution  to  the  integral  we  must  choose 


This  sets  a  lower  limit  to  K  lor  a  given 


a  value  of  E  smaller  than  E  . 

n  o 

u)  /d,  which  depends  on  the  altitude,  w  may  he  such  that  /j.-ii  is  smallest 
pe  1  r  r 

for  n=p  (say).  The  corresponding  lower  limit  for  k  is  fi  ed  by  reipiiring 

v  2  2 

that  E  <E  .  Next,  wc  note  that  .J  ,  -  .)  ,  is  an  oscillating  function  ot 

p  o  p+1  p-1 

the  argument  with  progressively  decreasing  amplitude.  Sin.  e  the  maximum 

contribution  of  G  ,  comes  from  the  vicinity  of  F~  I.  and  si  ice  G  ,  is  stronglv 
oh  o  oh  ■ 

peaked  at  E=E  ,  .1  positive  value  of  the  integral  (meaning  instability)  will  be 

0  l  n  11 

obtained  whenever  K,(l-E  \\  )  falls  within  the  positive  domain  of  J  ,  -  J 

x  p  o  1  p+1  p-1 

In  other  words,  the  regions  of  unstable  K,  will  alternate  with  those  of  stable 


The  largest  positive  value  of  the  integral  will,  of  c  nirso,  be  obtained 

1/2  1 
when  Kx(  1-1 /Eq)  is  equal  to  the  value  at  which  the  t  i  r  maximum  of  .l^+j  - 

d2  ,  occurs.  However,  for  the  same  values  of  u>  /it  and  K.  ,  which  ensure  the 
p-1  pc'  x 

positivencss  of  the  n=p  term,  there  will  be  some  or  all  11/p  for  which  the  value 

» 

J*"  ,  is  negative.  I'hesc  terms  in  the  sum 

may  cancel  the  positive  contribution  from  the  n=p  term  and  thus  eliminate  the 

instability.  However,  if  we  choose  k,.  in  such  a  way  that  li  ■  li  ,  but  E  NF  for 
7  »i  1  p  o  n  o 

all  n/p,  then  the  negative  contribution  from  the  n/p  terms  will  be  comparatively 
small  and  the  dominant  contribution  to  y  will  arise  from  the  n=p  term  in  the 


of  k  (1-E  /Fi  )  may  be  such  that  J*-  .  -  ./' 
xv  n  o  7  n+1 


sum.  Ill  is  sets  an  upper  limil  to  k^.  'flu-  opt  imal  value  of  k  ( ^  is  dot  ermi  lied  by 
the  width  of  G^flij.  Finally,  we  observe  that  the  upper  lubriil  mode  has  the 
largest  growth  rate  in  comparison  with  the  Bernstein  (cyclotron  harmonic)  modes. 


/ 


This  is  so  because  .'•»/.>«>  is  large  near  the  cyclotron  harmonics  (h  ,  n=2,3,4, 
...).  lienee,  we  shall  concentrate  on  the  upper  hybrid  instability.  In  this  case, 
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where  the  arguments  of  .1“  .  and  J“  ,  are  the  same  as  that  of  J‘  in  Eq.  (31). 
h  n+1  n-1  n 

In  order  to  obtain  values  of  y  we  now  consider  the  parameters  of  our 

5  -3 

interest.  In  the  lower  ionosphere  (100-170  km),  nQ  =  (0.6-3. 0)  *  10  cm  , 

B  =  0.35G;  then,  to  /ft  •  2.24  5.0.  Numerical  values  of  y  have  been  obtained 
o  Pe  . 

using  the  electron  energy  distributions  of  Jaspersc* .  The  collision  frequencies, 
v  ,  have  been  calculated  using  the  recent  cross  section  data,  and  the  model 
atmosphere  of  Jacchia11.  Electron  temperatures  were  taken  from  Jasperse1  calcu¬ 
lations.  The  collisional  damping  rates  are  given  by  Eq.  (25).  Results  of  these 
calculations  are  summarized  in  Table  1.  The  maximum  value  of  the  growth  rate, 

Ym,  corresponds  to  the  maximum  value  of  the  sum  appearing  in  Eq.  (32').  We 
notice  that  /  increases  with  altitude.  This  is  mainly  due  to  the  fact  that 
both  w  /ft  and  n(q1/n0  increase  with  altitude  (see  liq.  (32')].  The  integral  in 
Eq.  (52'|  decreases  with  altitude  due  to  broadening  of  the  distribution  function, 
but  the  decrement  is  very  small  over  the  altitude  regions  of  our  interest.  IVe 
have  verified  it  further  by  evaluating  the  integrals  with  G  j  =  6(E-Eo)/4tt  and 
then  comparing  these  with  the  actual  numerical  values.  The  collision  frequency 
and,  hence,  the  damping  rate  m  decreases  with  altitude.  This  is  due  to 
tion  of  neutral  species.  At  higher  altitudes,  however,  cl ectron-clectron  and 
electron-ion  collisions,  not  considered  in  this  work,  become  increasingly 
important . 

We  further  notice  that  the  value  of  ktv0/ft  for  which  the  growth  rate  is 

■}  ~>  i  >  •} 

maximum  increases  with  altitude.  .So,  b  k7  T  / m:2 “  =  (k"v7/)T  /2E  )  increases 

J-  e  X  o  e  o 

with  altitude  and  becomes  finite.  Consequently,  the  growth  rates  in  Table  1, 
which  are  obtained  for  b- *1 ,  have  to  be  modified,  from  the  dispersion  relation 
(19)  we  find  that  the  f i n i t e-b-mod i f icd  upper  hybrid  mode  has  a  frequency 
larger  than  its  previous  value.  l-'or  example,  when  u  /:>.  -  4.3  (140  km)  b~  1 . 0 
and  then  the  frequency  is  u>|j '5.4ft  whereas,  for  h^'l,  j*  4 .  1 4 1 ;  (see  Ref.  10  for 


the  solutions  of  the  dispersion  relation  when  I)  is  finite).  This  means  that 
instead  of  n=4,  n=5  term  will  contribute  most  to  the  integral  in  Eq.  (32'). 

The  growth  rate  will  be  smaller  than  the  one  reported  in  Table  1,  but  the 
reduction  is  small.  At  high  altitudes  electron-electron  collisions,  not  con¬ 
sidered  in  this  work,  will  be  more  effective  in  eliminating  the  instability. 

Finally,  at  altitudes  below  12S  km  the  collisional  damping  rate  exceeds 
the  maximum  growth  rate  of  the  upper  hybrid  mode,  and  so  the  instability  cannot 
be  excited.  As  we  show  in  the  next  section,  below  I21>  km  a  different  type  of 
instability,  namely,  the  electrostatic  electron  cyclotron  instability,  can  be 
excited . 


4.3.2  Electron  Cyclotron  Instability 

In  the  previous  section  we  considered  unstable  waves  with  small  but 
finite  k  .  Here  we  show  that  under  certain  conditions  waves  propagating 
strictly  across  the  magnetic  field  (k^=0)  can  also  he  driven  unstable  by  the 
suprathermal  electrons.  This  instability,  which  is  of  non  resonant  type,  is 
excited  at  the  electron  cyclotron  harmonics.  In  order  to  calculate  and  L)^ 
for  this  case,  we  change  over  to  the  spherical  coordinates  in  the  velocity 
space  and  after  some  straightforward  algebra  (see  Appendix i  obtain 
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In  deriving  F.qs.  (33)  and  (34)  we  have  put  k  -0.  The  complete  dispersion 


relation  is  then 
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n=l  m[(u+ivm)  -  n  \l  J  k^  y  hj 


where  N,  and  D,  are  given  by  F.qs.  (33)  and  (34),  a  (1  /•>)  exp(-b)  |  see 
h  h  n  n 

Eq.  (16)|,  and  we  have  denoted  k  by  kA .  Let  us  first  analvze  the  dispersion 

relation  in  the  limit  v  -0.  We  write  w=nst(l*x),  where  n-2  and  ■  1  (since 

m  — 

n  ,/n  1).  Then,  l.q .  (35)  becomes 

oh  o 


(4)  ) 


Then,  x  becomes 


X 


i 


(n2-l)ft 

nk  v 
x  o 


,2  2 
a  k.v 
n  J_  o 


nil 


(42) 


The  second  term  within  the  parenthesis  in  lit|.  (42)  represents  the  finite-b 
stabilizing  effect.  In  the  presence  of  collisions,  the  solution  of  the  complete 
dispersion  relation  (35)  is 


a)  =  nfi  +i 


Y  - 


(n“+l) u 


4n 


n>2 


(43) 


wiiere  y  =  nft(Imx)  and  the  additional  term  is  the  collision.il  damping  rate. 
The  maximum  value  of  the  collisionless  growth  rate,  denoted  bv  y  is 


m 


7  —.1/2 

a  k  \ 
n  m  ] 

a  / 


(44) 


Here  k  is  the  value  of  k.v  /U,  for  which  the  maximum  value  of  ll  is  obtained, 
m  J-  o  11 

Finally,  the  condition  for  exciting  the  electron  cyclotron  instability  in  the 
presence  of  electron-neutral  collisions  is 


(n4+l)v> 


'm 


4n~ 


(45) 


Although  accurate  determination  of  1  requires  numerical  integration,  it 
is  possible  to  obtain  an  analytic  expression  for  1  using  a  suitable  F  ^  that 
closely  represents  the  actual  distribution  function.  Such  an  analysis  is 
instructive,  if  not  so  accurate.  We  first  change  over  to  the  energy  representa¬ 
tion  so  that  Cq.  (37)  becomes 

I  .  4„Co  f  if  G  (E)  J2„P^  .  ,46) 

•o  *-  oJ 

We  then  choose 


l+« 


Goh(li)  =  4^T;  THTST  (r°  exP(-“L/F-o) 


(47) 
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which  satisfies  the  required  normalisation  condition:  4n  I  (H )  =  1. 

•o 

ft  can  be  easily  verified  that  such  a  distribution  is  peaked  at  E=E  ,  and 

®  > 

the  quantity  a  is  a  measure  of  the  width  of  the  distribution,  e.g.,  u=4 ( E^/A ) " , 
where  A  is  the  characteristic  half-width  of  the  distribution.  Our  model 
distribution  closely  resembles  the  ones  that  are  found  in  the  low  altitude 
regions  by  theoretical  calculation.1  Substituting  this  Gq^(F.)  into  Eq.  (46) 
we  find1'' 


(k  v  /ft)-  f(n+«) 
un'1r(l+a)r(2n+l) 


exp  (-k^v“/aft-)  M(n+1 -a; 2n+l ;k^v“/aft“) 


where  M  is  the  confluent  hypergeometr ic  function.  From  the  properties  of 

the  function  M  we  find  that  1  will  be  negative  when  (n+la)  <  0,  provided 
2  2  2 

(k^/v“/u<.’“)  •  1.  This  condition  tells  us  that  for  the  n-th  cyclotron  harmonic 

to  be  unstable  the  width  of  the  electron  distribution  has  to  be  such  that 

(A/2E  )  (n+1)  l/  i  .  e . ,  for  large  n,  the  distribution  function  has  to  be 

very  narrow;  otherwise,  the  mode  will  be  stable.  We  have  shown  earlier  (Eq. 

(41)1  that  the  value  of  unstable  n  is  determined  by  the  value  of  to  /ft.  Since 

pe 

'jo  /ft  increases  with  altitude,  requiring  larger  n  values  for  the  instability 
to  set  in,  and  since  the  widths  of  the  calculated  distribution  functions  also 
increase  with  altitude  (see  Fig.  1),  we  can  conclude  that  the  cyclotron 
instability  will  le  operative  only  in  the  low  altitude  regions.  The  approxi¬ 
mate  width  of  the  unstable  spectrum  is  given  by 


»“  (n  +  j)  / i 4a- 2) 


2 

(n  +  /  (4a-2) 


since  M  is  negat i  e  within  this  region.  There  are,  of  course,  other  regions 

of  negative  M  corresponding  to  higher  values  of  k  v  /ft,  hut  for  those  values  1 

will  be  .mailer  i  i  magnitude.  We  can  obtain  the  approximate  value  of  I  by 

taking  the  •*  •  •  1  I  init,  in  which  case  (i  ,  -•  A(E-E  )/4  n ,  and  liq.  (48)  reduces  to 

oh  o 

1  '  (5G 


Numerical  calculations  using  the  theoretical  distributions  of  Jaspcrse1  indicate 
that  Eq.  (SO)  is  in  excellent  estimate  for  1  in  the  low  altitude  regions  of 
our  interest.  For  illustration  purpose,  let  us  consider  the  n-4  mode.  If  we 
take  B  ^  -  O.isSC,  i  lie  value  of  id  /$»  given  by  Eq .  (41)  (with  2a^~l]  then 


corresponds  to  the  electron  density  ;it  an  altitude  ISO  kn. .  Irom  Fq.  (50), 
max|l|  =  max (-Jg) =0.23,  corresponding  to  k^v^/0  =  7.1.  Numerical  evaluation 
of  Eq.  (46)  for  n=4,  using  actual  G  ^(H),  yields  max|  1 1  =0.  18,  which  is  obtained 
for  k^v^/D  =  7.1.  The  agreement  is  even  better  at  lower  altitudes  where  (i  (E) 
is  narrower. 


In  so  far  as  the  stabilizing  term  in  Eq.  (44)  due  to  finite-b  effect  can 
be  neglected,  the  collisionless  growth  rate  of  the  electron  cyclotron  insta¬ 
bility  is  large  compared  to  that  of  the  upper  hybrid  instability.  In  the 

1/2 

former  case,  growth  rate  is  proportional  to  a  while  in  the  latter  case  it 
is  proportional  to  o,  where  o  =  n0j1/no‘^l.  In  Table  2  we  present  maximum 
collisionless  growth  rate  (without  the  finite-b  effect),  and  collisional 
damping  rate  as  functions  of  altitude.  The  altitude  is  dt  term  i  tied  from  the 
value  of  w  /S2,  which  in  turn  is  determined  by  the  eye  lot  ion  harmonic  n 
according  to  Eq.  (41). 

As  the  altitude  increases,  the  finite-b  stabilizing  effect  becomes 
increasingly  important.  This  is  due  to  the  fact  the  electron  temperature  as 

well  as  k  increase  with  altitude.  For  n=4  mode,  if  we  take  into  account  the 

m  '  1-1 
finite-b  effect  the  reduced  growth  rate  y  is  found  to  be  y  =11)  s  .  Above 

this  altitude,  the  stabilizing  term  dominates  and  so  the  modes  are  stable. 


4.4  SUMMARY  AND  DISCUSSION 

We  have  studied  the  electrostatic  instabilities  associated  with  the  E 
region  (below  170  km)  photoeleetron  distribution  functions  in  the  presence 
of  electron-neutral  collisions.  We  find  that  at  low  altitude  (at  or  below 
130  km)  waves  at  electron  cyclotron  harmonics  can  become  unstable  with  growth 

rates  larger  than  the  collisional  damping  rates.  As  altitude  increases,  this 

2  2 

electron  cyclotron  instability  is  extinguished  by  finite-b  (  k^T^/nurj  effect. 

At  altitudes  130  km  and  above,  on  the  other  hand,  upper  hybrid  instability 
can  be  excited  effectively.  Our  calculations  based  on  the  recent  theoretical 
electron  energy  profiles  of  Jaspersc*  show  that  the  upper  hybrid  instability 
can  operate  at  altitudes  as  low  as  120  km.  This  is  an  improvement  over  the 
calculation  of  Bloomberg4,  who  found  that  the  upper  hybrid  mode  is  stabilized 
by  collisions  at  or  below  160  km.  This  disagreement  witli  Bloomberg's  calcu¬ 
lation  is  primarily  due  to  the  fact  that  he  considered  the  collision  fre¬ 
quencies  of  l he  I  eV  electrons.  But,  as  we  have  shown  in  the  text.  1  he  damping 
is  mainly  due  to  the  Maxwellian  electrons  with  temperature  O.OI  cV  at  !.!<)  km. 
Consequently,  the  relevant  collision  I  r<  nem  y  is  an  null  a  of  magnitude 


smaller  than  that  of  Bloomberg.  Also,  our  calculated  collisionless  growth 
rates  are  somewhat  larger  than  those  due  to  Bloomberg.  This  can  be  attributed 
to  our  use  of  the  electron  energy  distribution  of  das perse1 ,  which  we  believe 
is  more  accurate.  At  high  altitudes  (above  170  km J  electron-electron  colli¬ 
sions,  not  taken  into  account  in  this  work,  are  expected  to  eliminate  the 
upper  hybrid  instability. 

Both  of  these  instabilities  arc  driven  by  the  suprathermal  electrons  with 
energy  peaked  at  4  eV.  Turbulence  arising  from  these  instabilities  is 
expected  to  produce  the  observed'*  • 3  anomalous  structure  in  the  electron 
distribution  function  in  the  2-4  eV  energy  range.  This  will  be  the  topic  of 
our  future  work. 
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5  APPROXIMATE;  ANALYTIC  SOLUTIONS  l-'OR  Tin;  AURORAL  PROTON  AND 
HYDROGEN  FLUXLS  AND  RELATED  QUANTITIES 

5.1  INTRODUCTION 

The  proton  aurora  results  from  the  precipitation  of  energetic  ( -  I  to  100  keVj 
protons,  with  the  peak  of  the  energy  spectrum  lying  below  10  keV,  into  the 
auroral  atmosphere.  Spectroscopically,  it  is  identified  by  the  readily  detect¬ 
able  Balmer  series  emissions  (H  ,  11  ,  and  11  j.  The  Balmer  emissions  are 

a  P  y 

radiated  by  the  moving  hydrogen  atoms  that  are  produced  in  the  charge-exchange 
collisions  of  the  protons  with  the  atmospheric  constituents.  The  proton  pre¬ 
cipitation  also  results  in  the  excitation  of  various  nitrogen  and  oxygen  emis¬ 
sions,  such  as  A3914,  A4709  N*,  and  A5577  01. 

Theoretical  calculations  of  energy  deposition  function,  ionization  rate, 

and  A3914  N_/I1_  intensity  ratio  were  initiated  bv  Chamberlain  | 19611.  Since 
l  p 

then,  Eather  and  his  co-workers  have  improved  upon  Chamberlain's  calculations 
utilizing  recent  experimental  and  theoretical  cross  sections  for  all  the  energy 
loss  processes  involved.  A  complete  list  of  these  works  can  be  found  in  the 
monograph  by  Val lance  Jones  [1974 1.  These  calculations,  however,  are  semi- 
empirical  in  nature,  without  the  detailed  knowledge  of  the  energy  distribution 
and  the  pitch  angle  distribution  of  the  precipitating  protons  as  a  function 
of  altitude.  As  Eather  [1967]  has  pointed  out  both  the  energy  spectrum  and 
the  pitch  angle  distribution  must  be  known  for  correct  theoretical  interpreta¬ 
tion  of  the  measured  hydrogen  line  profiles. 

In  this  paper,  we  calculate  for  the  first  time  the  energy  and  pitch  angle 
distribution  of  the  auroral  protons  as  a  function  of  altitude,  applying  the 
methods  of  linear  transport  theory.  The  proton  precipitation  problem  is  com¬ 
plicated  by  the  partial  neutralization  of  the  incident  protons  as  a  result  of 
the  charge-exchange  collisions.  This  leads  to  a  coupled  set  of  transport 
equations  -  one  for  the  proton  flux,  and  the  other  for  the  neutral  hydrogen 
atom  flux.  Another  related  complication  is  that  the  hydrogen  atom  resulting 
from  the  charge-exchange  collision  can  travel  large  distances  across  the  earth's 
magnetic  field  lines  before  being  converted  hack  to  a  proton  via  charge-stripping 
reaction.  Davidson  [1965]  used  a  Monte  Carlo  technique  to  analyze  this  trans¬ 
verse  diffusion  and  found  that,  for  5  to  20  keV  protons,  for  example,  an 
isotropic  incidint  pitch  angle  distribution  results  in  a  spreading  of  the 


precipitation  zone  over  an  area  as  wide  as  bOO  km.  Consequently,  the  relevant 
transport  problem  is,  strictly  speaking,  a  two-dimensional  one.  In  the  present 
work,  we  ignore  this  transverse  diffusion  and  obtain  approximate  analytic 
solution  ;  for  the  iroton  and  the  hydrogen  fluxes  in  a  plane-parallel  geometry. 

We  also  iresent  analytic  expressions  for  energy  deposition  rate,  and  ionization 
rate  as  a  function  of  altitude.  In  Sec.  ,  we  calculate  these  quantities  as 
well  as  he  densit  r  of  electrons,  which  result  from  ionization  and  the  strip¬ 
ping  of  il  atoms,  f >r  an  i sotropic-Maxwcl lian  incident  proton  flux,  using  the 
pseudoparticle  method  of  .Jasperse  and  Strickland  [197P]. 

5.2  CHARACTL-RISTI  1  FEATURES  OF  PRECIPITATING  PROTONS 

Fast  protons  -ntering  the  upper  atmosphere  undergo  charge- exchange 
collisions  in  which  the  incident  proton  captures  an  electron  and  becomes  a 
fast  hydrogen  atom.  The  H  atom  then  undergoes  stripping  producing  a  proton  and 
an  electron,  thus  completing  the  charge-changi ng  cycle.  These  processes  are: 

tl+  *  M  >•  11  +  M+  (1) 

11  -i  M  ►  ll+  -i  M  +  e"  (2) 

where  M  denotes  any  atmospheric  constituent.  The  cross  sections  for  these 

-16  2 

processes  are  quite  large  (>10  cm‘) .  As  a  consequence,  an  initially  pure 
H+  flux  ;oon  becomes  a  mixture  of  ll+  and  11,  and  after  only  a  few  charge- 
changing  cycles  ar  equilibrium  flux  is  established  in  which  the  fractional 
composition  is  a  function  of  energy.  For  auroral  energies,  this  equi 1 ibriation 
occurs  well  above  the  altitude  where  significant  energy  losses  due  to  ioniza¬ 
tion  anil  excitatiin  of  the  atmospheric  constituents  can  occur. 

In  addition  to  processes  (1)  and  (2),  charge-changing  processes  such  as 
II*  <■  M  >  if  •  M* +  (3) 

If  »  M  >  II*  ■  M  +  2e  (4) 

can  occur.  However,  the  cross  sections  for  these  processes  are  very  small  and 
the  li  content  of  the  equilibrium  flux  is  'I"  in  the  auroral  energy  range, 
because  lonizatioi  and  excitation  processes  involving  11  have  cross  sections 
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comparable  to  those  of  H  and  II  processes,  it  is  a  pood  approximation  to 
neglect  H  in  the  proton  auroral  analysis. 


5.3  BASIC  EQUATIONS 


For  a  single  constituent  atmosphere,  the  linear  transport  equations  for 
energetic  protons  and  hydrogen  atoms  in  a  plane-parallel  g.-ometry  are 


^3“  +  n(z)|op(EJ  +  o®x(H)  +  o^*(E)  +  o 1 ° (E)  j.  j 4^ (z , I  . p 


=  27in(z)  £  |dE^d^-£  o  .  ►  F ,  u  J  <J>  (z,L',m") 

k  p.i  p 


:/ 

ij* 


+  2in(z)  T.  I  dli'dy '  a'j’Vi'.u'  *  E,p)  *„(z,i:',y 

:  I  . 


(5) 


[u  37  +  n(z)|oJ(E)  ♦  o®X(E)  +  ajj£(0)  +  o01(E)IJ*jj 


(Z,d,p) 


=  2irn(z)  £  |dE-dn"  £  0^(1;', p'  •>  E,p)  4>H(z,E'  ,p') 


f 

if 


10, 


+  2pn(z)  Z  J  dE  'dp  '  a.  (E',p'  -*•  E,p)  $  (z  ,  E ' ,  p  ') 


(6) 


where  <t>'s  are  the  particle  fluxes,  z  is  the  altitude,  p  is  cosine  of  the  angle 

between  the  particle  velocity  and  the  positive  z-axis,  whi  h  is  parallel  to 

the  geomagnetic  field  lines,  and  n(-)  is  the  neutral  density.  Here, 
k 

o.(E',p'  -*  E,p)  is  the  differential  cross  section  per  unit  range  (referred  to 

hereafter  as  just  the  differential  cross  section)  for  collision  between  the 

precipitating  particles  and  the  neutral  particles  in  which  the  neutral  particle 

makes  a  transition  from  the  ground  state  to  the  final  state  j.  The  summation 

index  k  labels  elastic  (ef ) ,  excitation  (ex),  and  ionization  (i)-type  collis- 

k 

sions.  The  corresponding  total  cross  sections,  denoted  In  o  (E) ,  is  related 
to  the  differential  cross  sections  by  the  formula 


ok(E) 


JdE'dp 


a.  (E',p' 


E  ,p) 


C) 
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The  subscripts  P  and  11  represent  cross  sections  associated  with  the  protons 
and  the  hydrogen  atoms  respectively.  In  addition,  is  the  differential 

cross  section  for  the  charge-exchange  process  (1)  in  which  the  ionised  neutral 
particle  is  left  in  the  state  j,  and  o is  the  equivalent  differential  cross 
section  for  the  stripping  process  (2).  The  corresponding  total  cross  sections, 
u10(E)  and  o^(E),  are  related  to  the  differential  cross  sections  by  the  formulae 
simi lar  to  Eq .  (7) . 


In  writing  Eqs.  (5)  and  (.6)  we  have  ignored  any  electric  field  that  may 
be  present  in  the  auroral  atmosphere  and  assumed  that  the  geomagnetic  field  is 
uniform.  Under  these  assumptions,  protons  lose  energy  only  via  collisions  and 
such  collisions  result  in  discrete  energy  losses. 


5.4  FORWARD -SCATTERING'  AND  AVERAGE  DISCRETE  ENERGY-LOSS  APPROXIMATIONS 

The  differential  cross  sections  for  all  the  processes  involving  protons 
and  H  atoms  are  very  strongly  peaked  in  the  forward  direction  for  incident 
energies  above  a  few  hundred  eV  [McNeai  and  Bircly,  1975].  This  suggests  that 
except  toward  the  end  of  the  precipitation  one  can  make  the  forward  scattering 
approximation.  We  shall  use  this  simplifying  approximation  with  the  under¬ 
standing  that  our  analysis  will  be  only  valid  for  E>E  .  . 

J  min 

The  loss  function  M.E) ,  the  energy  loss  per  unit  path  length  per  molecule, 
for  a  particular  inelastic  process  (excitation,  ionization,  charge-exchange 
or  stripping)  may  simply  bo  defined  as  the  product  of  the  average  energy  loss 
W  and  the  total  cross  section  associated  with  the  process.  Edgar  et_  aj^  [1973] 
have  given  the  loss  function  for  each  inelastic  process  as  a  function  for  each 
inelastic  process  as  a  function  of  proton  energy.  Using  the  cross  sections 
given  by  McNeai  and  Bircly  [1973],  and  assuming  equal  average  energy  loss  for 
both  protons  and  Ii  atoms,  wc  have  calculated  W  as  a  function  of  E  from  the  loss 
functions  of  Edgar  ot_  tpl_  [1973]  and  find  that  it  is  a  weak  function  of  energy 
in  the  auroral  energy  range  and  that  it  is  almost  the  same  for  all  the  inelastic 
processes.  In  the  average  discrete  energy-loss  approximation,  IV  is  taken  to 
he  a  constant  (=  28  oV)  and  the  same  for  all  the  inelastic  processes. 


With  these  approximations,  wc  write 

nc5'(i:\ti'  •  E,n)  -  (2ti)_1  oC?'  (E ' )  S(E'-E)  <Un'-|,)  (8) 
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E  aeX(lT,M'  -  E,p)  =  (2n)~J  E  u?x(E')  A  f  I: (E+Kj  |  «'.(u  -n) 

j  J  J  J 

e  01(f',p'  -  i:,pj  =  (2m)'1  >:  o!(ir)  a  L  r: (i..wi  |  .s  (u 
j  J  j  J 

E  ck°(K  '  ,\i  '  *i:,p)  =  (2n)  1  >:  oJt)(l:')  All- *-(li+K)  |  .S(|i'-.J 

j  '  j  J 

E  01  (E n '  -  E,U)  =  (2ir)_1  E  o01(E')  «5[E'-(E+W)1  6(p '-p) 

j  J  j  J 

for  both  protons  and  H  atoms.  Here,  E  oPx(E')  is  the  total  excitation 

j  J 

0X 

section  a  (E'),  and  so  on.  Substituting  Eqs.  (8)  through  (12)  into  Eq 
and  (6)  we  obtain 

[p  ±  +  n(z)  0p(E)J  *p(z,E,y) 

=  n (z)  |apX(E+W)  +  op(E+W)j  <fcp(z,E+W,p) 

+  n(z)  o01(E+W)  <t>H(z,E+W,p)  , 

and , 

[m  ~  *  n(z)  oh(E)]  <t>u(z,E,p) 

=  n(z)  |^o®X(C+W)  +  oJCE+HOj  *n(z,[->tf,y) 

+  n(z)  o^(E+W)  ^(ZjE+W.p) 

Here  we  have  defined 

a  (E)  =  aeX(E)  +  o1  (E)  +  a10(E) 

P  P  P 

nH(E)  5  o**(E)  >  O^E)  ♦  o01(E) 


:ss 


(9) 

(1«>) 

(H) 

(12) 

cross 
s.  (5) 

(13) 

(14) 

(15) 

(lb) 


and 


Wo  not  id1  lh.it  the  clastic  scattering  terms  drop  out  as  a  result  ol  the 
forward-scatter i ng  assumption,  and  tlie  coupling  between  the  proton  ami  II 

atom  fluxes  arise  due  to  the  charge-changing  processes  represented  by  o  and 

01 

a 

5.5  SOLUTIONS  FOR  THE  FLUXES 

Since  the  equi li briat ion  of  the  incident  proton  flux  occurs  in  the  altitude 
region  where  energy  loss  is  unimportant  due  to  very  low  atmospheric  density, 
we  can  analyze  the  prec i pi  tat  ion  problem  in  two  stages.  First,  we  neglect  all 
energy  losses  and  show  that  equi 1 ibration  due  to  charge-changing  processes  is 
a  direct  consequence  of  Lqs.  (13)  and  (14).  V.'e  calculate  the  fractional  com¬ 
positions  of  protons  and  H  atoms  in  the  equilibrium  mixture.  Then,  we  solve 
t Ire  full  transport  equations  for  the  fluxes  taking  the  equilibrium  fluxes  as 
the  given  boundary  values. 

5.5.1  I  qiti  1  i  hr iat  i ng  Fluxes 

If  we  neglect  all  energy  losses,  Eqs.  (13)  and  (14)  reduce  to 


[u  •"  +  »U)  'J10(K)j  Tp(z,l:,n)  =  n(z)  o0I(K)  4>H(z,E,u)  07) 

n  -tr  +  n(z)  uOI(l:.)J  4([(z,l.,u)  =  n(z)  «J°(E)  4>p(z,E,u)  (18) 

Adding,  the  two  equations  we  obtain 

--  |t  (z.E.u)  4  !.,(-, E, u)  1=0  ,  (19) 

p  II 

whi  cli  i  ittp  ties 

■t'pU.n.n)  +  i’(,f z,i:,ii)  -  (-0) 


where  ‘l>()(E,n)  (for  -1  p •  0 )  is  the  proton  flux  incident  (at  z=<”)  on  the  atmo¬ 
sphere.  liquation  (20)  is  just  u  statement  of  the  conservation  of  flux. 

Inserting  liq.  (2l>)  into  either  ( 1 71>  or  (18),  and  defining  an  equ  i  1  i  br  i  at  ion 
depth  r  as 

d  -  -  (  '  **'  *  )  n  (7.  )  d z  (21 ) 


3(> 


I01 

‘i’pfTp  (z  ,E)  ,  1: ,  m)  =  ^qCEjM)  jq  qj" 

a  +  a 


1  + 


10 

1 

To 


t'xl>(  t,/f) 


(22) 


10  ,  x 

,E,P)  =  4>0(E,m)  -Jq^ - qT  (l  -  exp(i  /p)J 

n  o  +  a  '  / 


(23) 


for  -l<u<0,  and  zero  otherwise,  since  ( l. ,  ii )  =  0  for  (t-pl,  i.e.,  no  particles 
are  incident  from  below.  According  to  F.qs.  (22)  and  (23),  the  fluxes  approach 
their  respective  equilibrium  values,  given  by 


01 

4>ec^ CE  u)  =  _ ^ „ 

P  10  01  0 

a  +  a 


4>_ (R , p)  ,  for  -lp<0  , 


(21) 


and. 


10 

^(E.vO  =  10  oi  To 


$  (E,u) ,  for  -l<(i<0  , 


(25) 


a  +o 

asymptotically. 

For  an  isotropic  incident  proton  flux,  the  hemi spher i cal ly  averaged  values 
of  the  equilibriating  fluxes  are 

r 0 

4>p(Te(z,E)  ,E)  -  2tt  J  dp  <5>p(Tg(z,E)  ,E,p)/2n 


=  *pq(E) 


h 


10 

01  ti2^Te('Z,L^ 


(26) 


and , 


/. 


0 


4>H(Tetz»Ii)  T)  -  2n  J  dy  ‘?>ll('l0  (Z»E)  ,E,u)/2tt 


-1 

^  *eclrr 


H  e 


=  (E)  [1  *  r2(Te(z,E))] 


(2?) 


Here, 


ipq(E)  f  du  $pq(E,u)  , 


(28) 
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(29) 


*[,q(ii)  /  dp  , 


and,  L:,(x)  is  the  second  exponential  integral  function, 
sphere,  the  most  effective  charge-changing  process  are 


In  the  upper  atmo- 


H+  +  0  •>  H  +  0+ 


,  10, 
(a  J 


H  +  0  +  H+  +  O  +  c'  (o01) 


due  to  large  abundance  of  oxygen  atoms.  We  have  calculated  values  of  4>p  and 
as  a  function  of  altitude  for  an  isotropic-Maxwellian  incident  proton  flux. 


■I'n (Ci.)  *  ' ’-T\  I 

V- 


— li  ex|i(-l:./ho)  ,  for  -l‘_y<0  , 


where  C>  is  the  total  energy  flux  in  the  downward  direction  and  E  is  the  mean 

^  1001  ° 
proton  energy,  using  o  (E)  and  a  (E)  of  McNeal  and  Birely  [1973]  and  using 

a  model  atmosphere  [Jacchia,  1977;  1000°K] .  We  find  that  for  Eq  =  8  keV,  for 

example,  the  fluxes  attain  their  equilibrium  values  at  approximately  290  km, 

i  ,c. ,  Z  ~  290  km. 
eq 

In  the  following  precipitation  analysis,  we  shall  take  z=zCq  as  the 
boundary,  and  solv  ■  f-qs.  (13)  and  (14)  wi  th  <t>p^  and  ,  given  by  Eqs.  (24) 
and  (23)  respective  ly,  as  the  fluxes  incident  on  the  boundary. 

5.3.2  Precipitating  Fluxes 

Trans  forming  '.qs.  (13)  and  (14)  into  equations  in  terms  of  the  'optical' 
depth  . ,  where 


/“eq 

n(z')  dz'  , 


we  obtain 


f  3  V  O  (n+W) 

V*  ~  -  lj  ♦  p(t(z,E),E,jj)  =  -  0^(1;)  ♦,.(•■(=,  r*W),l>W,,.) 

-  ^H(r(z,r.,W)  ,E+W,p)  ,  (32) 


and, 


f  a  \  o  (E+W) 

W  37  -  lj  ♦jjfTtz.Ej.E.u)  =  -  ~jny  ♦|1(T(z,E+W),l->Wtll) 


0  oiw  ♦P(t(z>I->w),F-*W,p)  .  (33) 


Here,  we  have  introduced  op  £  op  +  up  ,  Oj (  +  for  brevity  of  notations, 

and  assumed  op(E)  =  Oj((E),  which  is  a  fairly  good  approximation  over  the 
auroral  energy  range. 


The  formal  solutions  of  Eqs.  (32)  and  (33)  are 


*p(t,E,u) 


<t>pq(E,p)  exp(t/u) 

V°+W)  fT  dt 


/ 


opCE)  J  u  Vbnt,E=W’M)  oxPf  (r-t)/M] 


/  ^^H(but,E+W,M)  oxp[(r-t)/pl  , 


and 


4>h(t,E,u)  =  4>}eq(E,p)  exp(T/p) 
oH(E+W)  .t  , 

/■  -  oprE)  Jo  r$H(V’c+w>lJ)  ‘'xpf(i-n/lJl 
'  gup((i~)W)  /o  r* Pcv«E+w*^  exp  t CT-tj/uj 


(34) 


(35) 
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for  -hu<0.  Here  b(|(l:.)  Op(l;+W)/o  (E).  Adding  Eqs.  (34)  and  (35)  we  obtain 
4'p0  +  »1I(T,l:,p)  =  4>q(1;-,m)  exp(x/p) 

-  I»U(K)  J  ~  pp(but,E*W,u)  *  *H(but,E+K,w)J 

x  exp  l  (x-t)/u|  (36) 

under  the  assumption  °p=0|j-  I11  obtaining  Eq.  (36)  we  have  used  Eqs.  (24)  and 

(25).  This  equation  for  the  combined  flux  can  be  solved  by  iteration  procedure, 
which  is  equivalent  to  the  multiple-scattering  method  used  by  Jasperse  and 
Strickland  (1979)  in  solving  the  auroral  electron  precipitation  problem,  and 
the  complete  solution  is 

(O 

<t>p(r>E,u)  +  *1((t,1;,u)  =  •l'0(E,u)  exp(T/u)  +  ?  bnJ(E) 

n=l 

x  .D0(!>nW,u)  Jh^t/p.E)  -  UJO.E)  exp(r/w)J  ,  (37) 

for  -l'u  0,  and  zero  otherwise.  In  the  above  equation, 

HjfX.E)  =  au(l.)  oxp(l>ux)  ,  (38) 

and  the  higher  order  11  functions  can  be  obtained  from  the  recursion  relation 

exp  I  - x)  ll|i+1(K,in  -  Hn+1(0,i;)  = 


dx  •  exp(-x  )  |lln(l.  x  M:+W)  -  If  (O,l-:+l0  exp(h^x ' )  | 

so  that 

H.IX.E)  -  aJ.a,1  exptb^x)  -  cxp(b  x)  , 

H  IX.ii)  =  a  a  a...  oxp(l>  xj  -  a.  a  a,.  exp(l>  x) 

-•»  .■>.*>  .*>  1  1  3 i  v*).*)  .*>2  11  1  u 

'  ;,.,;ia72;,21  *  a33a22'1 1  I  t-xl,fl\xx>  ’ 


(39) 


(40) 


(41) 
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etc.  In  the  above  formulae  wc  have  introduced  the  following  definitions 


bmn  =  bmn(E)  "  op(E+mW)/0p(E+(n-l)W)  , 


a  =  a  (E)  =  a  (E)  s  fl  -  b  (E)l  . 
mn  mn  mn  L  run  -> 

That  Eq.  (37)  solves  (3(>)  exactly  can  be  verified  by  direct  substitution  and 
using  the  recursion  relation  (39). 

Substituting  <|>  from  Eq .  (37)  into  Eq.  (31)  and  doing  integration.,  we 
obtain 


^p(t,E,a)  =  <J>pM  (E  ,  vi)  exp(  i/p) 


+ 


a01 (E+W) 
Op (E+W) 


<JO 

2  <J>o^,'+nW’11  *bnl  ( *  )  [fln(  ‘/l1  .E)  ll(i  (0,  E  |e\|i(  i  /  n)] 
n=l 


op (E+W) -n01 (E+W) 
Op (E+W) 


b  ( !  •:) 

It 


x  exp[(T-t)/uj 


MD 


In  doing  the  integrations  we  have  used  the  recursion  relation  (39)  and  the 
definitions  (42)  and  (43). 


We  now  solve  Eq.  (44)  by  iteration,  where  t  lie  cerotli  order  iterate  is 

oo 

o  1 

.  iO  exp(i/y)  +  2*  4-()(l:.+nK,n) 

'  '  n= 1 

X  hnl^  [Hn(T/M,E)  -  lln(0,E)  oxp( t /g )]  (IS) 

After  successive  iteration,  the  resulting  solution  can  be  expressed  in  a  closed 
form  as: 
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:r>' 


win  ri¬ 


ll:  I 

I 

so  t  hat  , 

(K> 


HO  y 


■ )  -  t-lo)  I  i 
I* 


!  tlcr+nW,iO  !)nl  (10 


■h' 


(t/n 


♦  2l  cr)  U! 

nl  "  LP 

ou 

01 

,10-11  (0,10  exp(i/u)l  +  - - [■ :+  X  4'  (E+mlV, 

"  1  a  p(L+  (n+1  )W)  oK 


(ot 


X  lm[(l':){,1m(r/lJ’L)'Hn,{0’E)  exPfT/^} 

is  g i vi  n  by  Fq .  ( IS)  ,  anti 

>1 


v) 


(46) 


(1.+1- ) -o01  (L+W) 


L 


o  (F  +  K) 


J 


o  (L+2W) -cs01  (F+2W) 

o  (H+nW)  -a01  (F+nlV) 

o  (F+2W) 

P  J 

a  (E+nW) 

_  P  J 

(4- 


11 


Cr+w  i  ->° 1  (i>wi 

0  IF+IV) 
i1 


o  li-:+v  )->i01  (i:.+w) 


a  (E+W) 


a  (F.+  2W)  -  o01  (F.+2W) 


“o  (1:  +  2W) 

P 


OtC  . 


Fciuations  (57),  (41)  and  (47)  constitute  the  complete  solution  for  the  preci¬ 
pitating  proton  ant  hydrogen  fluxes  subject  to  the  approximations  that  were  made 
in  the  course  of  tie  analysis. 

A  few  comment-  about  the  solutions  are  in  order.  The  linear  dependence 

of  the  fluxes  at  u: I  altitudes  on  the  incident  flux  is  evident.  The  solutions 

yield  no  hackscatti  red  flux,  i.e.,  4>=o  for  0^ p <  1 .  The  reason  for  this  is  that 

in  the  forward-scattering  approximation  elastic  scattering  terms  drop  out  and 

no  particles  are  scattered  into  the  backward  direction,  bach  factor  in  (R)  is 

n 

a  very  small  quantity  (-0.00)  over  the  auroral  energy  range,  so  that  the  sum 
in  Fq.  (4(>)  convenes  estremely  rapidly.  For  practical  calculations,  only  the 
n=l  term  in  the  sun  should  he  quite  sufficient. 
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5 . 6  HEM1SPHERICALLY  AYERAliHD  Rl.lJXl:.S,  INI.IU'.Y  DEPOSITION  RATH 
AND  IONIZATION  RATI! 

In  this  section  we  give  analytic  expressions  for  various  measurable 
quantities  for  an  i  sot  rop  ic -Maxwe  I  1  i  an  incident  proton  flux  as  given  in  hq.  (30) 


5.6.1  Hemispherical  ly  Averaged  I- 1  uxc 

,u) 


$  (  '•  ,E)  I  dp  4>_  C  r » I'- .  1 

p  i\  p 


n  1 


a01  (l.+nW ) 


o^(H+nW)  +  ''()1(|>nW) 


( H  +  nlV) 


*|Kn(t,i:j  -  11^(0,11  I. ,( i  )  j 

•ft)  %  -4 

x  exp(-  (K+nW)/liQ)  b^K  ^t  ,11)  -  IMO.E)  I  , (  1  ) | 
01  _°°  _ 


m=n+l 


(18) 


and,  according  to  liq.  (37), 


Vt.F.) 


^Eq 


j  E  exp (-Ii/Eg)  E9(t) 


fXJ 

•I 


2,(E+nW)  exp(- (E+nW)/F.  )  h  j  <K  (t.E)-II  lO.li)  I;-2 ( i  )}• 
n=l  ‘  "  - 


-  ^(t.E)  , 

where  E2(x)  is  the  second  exponential  integral  function  defined  by 


(49) 


i 


43 


MX) 


Jt  t  exp(-xt) 


(50) 


and , 


f° 

KnCX,U3  J  dp  IIn(x/p,F.) 


(51) 


It  can  be  easily  verified  that  the  K- functions  are  the  same  as  the  11-functions 
with  the  exponential  functions  replaced  by  the  E^-functi ons  everywhere. 


5. (>.2  Energy  Deposition  Kate 

The  energy  deposition  rate  n^fz),  defined  as  the  energy  deposited  per  unit 
volume  per  unit  time  by  the  precipitating  fluxes,  is  given  in  plane-parallel 
g comet  ry  by 


Oj;  (  7  ) 


2  n  n  ( z  : 


f  *  f 

Jo  -1 


dp  p  E  0  (E)— 

p  t)  T 


l<J>p(T,E,p)  + 


4,|j(t,F.,p)  ] 


(52) 


recognizing  that  z)/.')z  n  (z  )o^  (E)  D4»  f  i  )/3t  .  With  the  aid  of  Eq.  (37)  this 
becomes 


n,.(r)  =  ( 


(t) 


(-f)  '•  (  z  )  f  dE  E  o  (1.)  E  exp (-1I/E  .)  E, 

i:;  e  .  *  L  ' 

0  in  i  n 


+  (Em.W)  exp(-  (E+nW)  / 1!  Q )  '^{Vt.HI-H/O.E)  (t  )| 


(55) 


where 


l.(X.li)  /  dp  p  ---  II  1 x/p ,E ) 

1 1  ^  |  ^  A  Jl 


I) 


a 


(54) 


The  rule  for  1,  (X,E)  is:  find  'll^fX  ,E)/5x  and  then  replace  t  he  exponential 
functions  by  the  E  ,-fuuct i oils  everywhere.  As  discussed  in  the  beginning  of 
section  5.4,  our  analytic  solutions  for  the  fluxes,  based  on  the  forward- 
scat  ter  i  ip*  approximation,  are  only  valid  for  li>l;  .  .  In  doing  the  energy  integral 
in  Eq.  (55)  from  I.  .  to  ,  instead  of  from  0  to  ■»  which  defines  n,.,  we  lose 
contribution  to  the  integral  amounting  a  few  percent  for  the  auroral  protons. 
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5.6.3  Ionization  Rate 


An  important  quantity  describing  t  he*  proton-hydrogen  energy  deposition  in 
the  auroral  atmosphere  is  the  ionization  rate  rr(z),  defined  as  the  number  of 
electron-ion  pairs  created  per  unit  volume  per  unit  time.  In  proton  auroras, 
in  addition  to  the  impact  ionizations  by  fast  protons  and  hydrogen  atoms,  an 
electron- ion  pair  is  created  at  the  end  of  each  charge-changing  cycle. 


In  plane-parallel  geometry,  n.(z)  can  be  calculated  from 


\(z)  - 2” "(!)  f  j 

^min  * 


dq  u  P.  a  (E)  —•  1 4>  (T,l',uV<t>[j(T,i;,u) 


(55) 


where  Wg(E)  is  the  average  energy  expended  in  creating  an  electron- i on  pair  due 
to  impact  ionizations  and  charge-changing  processes.  Inserting  F.q.  (37)  into 
(55)  we  obtain 

h  0p(E)  ^E  exp(-E/E0)  L;2(t) 

2(E+nW)exp(-(E+nW)/E0)  bnl  {1^(1  ,E) -Hn (0,0) 02 (i )  jl  .  (5<0 


The  underlying  assumption  in  writing  Lq .  (55)  is  that  the  energy  deposited  by 
the  precipitating  proton-hydrogen  fluxes  is  entirely  expended  in  producing 
electron-ion  pairs.  This  is  a  reasonable  assumption  for  I!  E  .  1  kcV,  since 

the  energy  loss  due  to  excitation  processes  is  comparatively  small  owing  to 
small  cross  sections. 


We  take  for  Wg(E)  the  values  calculated  by  Edgar  et  a  1  [1973].  Typical 
values  are  30,  27,  26,  27,  28  and  32  in  eV  for  1,  2,  5,  10,  20  and  50  keV  proton 
energies,  respectively.  Because  of  charge-exchange  a  great  number  of  ions  are 
created  at  low  energies  leading  to  a  dip  in  W.(E)  below  10  keV  preton  energy. 

We  further  point  out  that  in  evaluating  the  energy  integrals  in  Eqs.  (53) 

and  (56),  the  integrals  may  be  truncated  at  some  E  .  When  this  is  done  the 
°  1  max 

infinite  sums  terminate  at  a  maximum  value  N  given  by  N  =  (E  -E  .  )/W. 

*  max  mm 
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5.7  PSHUDGPARIICLI.  MITl.OD 


In  the  formulae  of  si-cticiis  5.6.1  through  5.0.3,  VV  is  50  cV  and  the 
incident  auroral  proton  energies  range  from  1  to  00  keV.  This  means,  the  sums 
in  the  formulae  extend  to  as  many  as  2,000  terms  representing  as  many  as  2,000 
scatterings  before  an  energetic  proton  loses  all  its  energy.  However,  a  good 
approximation  to  these  large  sums  may  be  obtained  by  introducing  the  notion  of 
pseuuoparticles  |J: asperse  and  Strickland,  1970],  A  pseudoparticle  is  a  particle 
which  has  a  cross  section  W/W  times  smaller  than  the  real  particle  but  has  an 
average  energy  loss-  per  inelastic  collision  V/ IV  times  greater,  where  W>>W.  In 
describing  the  energy  precipitation  in  terms  of  the  pseuuoparticles,  we  simply 
replace  o  by  o  and  W  by  W  everywhere  in  Has.  (48)  through  (56),  i.e. 


o(r.)  -  o  (E)  (Vl/W)  o  (E) 

(57) 

W  >  w  (W/W)  w 

(58) 

lie  notice  that  under  such  transformations  the  energy  loss  function  L(E)  Wo(E) 
remains  unchanged;  hut,  the  number  of  terms  in  the  sums  is  reduced  by  the  factor 
W/K,  i.e.  N  VN (W/W) 

Using  the  pseudopart ic le  method  we  have  calculated  the  various  quantities 

in  sections  5.6.1  through  5. (>.3.  We  have  found  that  a  converging  answer  can 

be  oV.tr Lned  with  a  suprisingly  small  number  of  pseudoscatterings.  For  example, 

tor  an  i  sot  rrpic-Maxwe  1 1  i  an  incident  proton  flux  with  F.  =  8  keV  only  10  pseudo- 

scatterings  (i.e.  U-IO)  spanning  the  energy  range  from  H  ^  =  2.5  keV  to 

H  -  57.5  keV  yield  good  answers, 
ma  x  b 
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We  change  over  to  spherical  coordinates  in  velocity  '.pace  using 

» 

v  =  v  cos",  vA  vs  imj,  v  dv^dv,,  v“dv  sino  do 

31-  ,  v  al-  .  31-  v,  3I; 

oh  _  >•  oh  .  _ <vh  oh 

3vn  v  3v  ’  3va  v  3V 
Then  the  expression  for  becomes 


Nh  = 


2Tte(io+i 


I  /■'  -  *  r 


d()  s  i  IV 


J“  Ik  vs  i  110/S.’) 
_n_ _i_ 

s,+  i\>  -k,.vcosf!-n:j 

h  o 


31 


dv  V 


3v 


taking  into  account  the 
With  the  substitution 


identity  X  »J“ 
n  n 


l.et  us  consul  r  the  first 


1 


ui+iv.  -k.,vcos0-nS/ 

h  a 


I 


dr  exp  I  - i ( m+ i v 


ll-k^VCOS". 


in: )  i 


and  using  the  addition  theorem 


+  On 

2 

J  (a  sin  b)  =  ,  .1“  (a/2 )  cxp(2inh) 

o  *— <•  n 

we  can  carry  out  the  integration  with  respect  to  o .  The  result  i: 

^ .  f  oh  /  ,  sin  Av  r  •  /  •  v  i 

I  =  -2i  I  dv  —  —  J  di - —  oxp{  -  i  (oi+  I  )  i  I 

Ju  ‘  -™> 


2  t  ?  2 

where  A  =  [k^)  +  ( 4^/11“)  sin  Sit/21  .  Integrating  by  parts  with 

to  v  wc  obtain 

I  =  2  i  j  dv  r  )h  ( v )  /  dr  cosAv  exp{ -i (w+i h  I 


(Al  ) 


integral . 


( A2  I 


respect 


i  Ad  1 


4!) 


In  the  case  k#.  =0 ,  we  can  use  the  identity 


V 

cos(  (2k^Vy/S2)  siiUCi/.M  =  ^  J.,  (ik^v/ii)  exp(inS2i) 

n  =  -“' 

and  integrate  with  respect  to  t.  The  result  is 


(A4) 


1  =  -2 


2  J 

»» —  •',1 


J2n(“V/u) 

dv  !■  ,  (v)  - 

_  oli  (i'+ 1  \>.  -mi 

n=-"'  ■'o  h 


(A5J 


Using  (A.ri)  in  (  Al)  we  get 


N 


•‘-"-e  f  dv  !•'  ,  (V) 
U  HI  J  oil 

A) 


;  ,  -  ,  V 

1  -  ll..+  l^(1J 


n  -- 


j  \\  -nS2 

h 


Pol  lowing  l  he  procedure  leading  up  to  laj.  (  A2J  ,  wo  <;ct 


(AO) 


J  ii  \) 


f 


n  J  dv  v  |;0h(v)  /  dt  k'X|>l-i  Ov-i  v,Jtj  CAT.) 

Oil  ■'ll 


V 


rv 

Writing  sin  Av/A  -  I  dv  eosAv  and  performing  the  integration  with  respect 

■'ll 

to  i  using  liq.  (Al  )  ,  we  obtain 


li,  I 


f..,  st 

Z  1  ,jv  y  \:  (-yj  f  ll V  .1  (2k  v'/i!)  (AS) 

n'+ 1  v.  -m2  oh  I  2 n  jl 

oil  n  =  -“>  li  Jo  Jo 


(•'ina  I  I y ,  us  i  ug 


f  dx  '  .1  ( x  ')  -  .1  , 

Jn  -n  Z  2n+Jl,+  l 

P=‘> 


we  obta i n 


(x) 


II V 


k  n 


V 


s,  oh  n 


(u+  i  v,  -  ni2 


f  dv  v  I  ^  J  (2k  v 

J  oh  i '  2n+2pi  I  i- 

p- o  1 


,/d)  (A1.)) 


r>d 


TAlii.i:  1 

Maximui:i  collisionless  growth  rate  y  ,  collision.il  damping 

rate  w.,  and  other  parameters  of  the  upper  hybrid  instability 

as  a  function  of  altitude,  for  B  -  0.35(’i.  Here,  k',  (  k.v  /:  ) 

o  o 

and  KA/K1|(=k. /k,,)  refer  to  the  values  for  which  maximum  growth 
rate  is  obtained. 


120  3.22  2.08  x  10  1.5  x  10'  8.  )<>  x  10“  1.10  x  10  8.8  ,  10.5 


Altitude 

(kmj 

u>  /Si 
pe 

120 

3.22 

130 

3.68 

140 

4.3 

v  (s 

! 


m .  ( s  ) 


Y  (s  ) 
m 


K  K  /k 
x  i  x  n 


tabu:  2 

Maximum  col  1 i sion less  growth  rate  Ym>  collisional 
damping  rate,  and  other  parameters  of  the  electron 
cyclotron  instability  as  a  function  of  altitude, 


for  B  =  0. 351! . 


